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ABSTRACT 


The program described in this report enables a user to 
compare the values of meterological parameters derived 
from data obtained by the SMMR instrument on NIMBUS-7 with 
surface observations made over the ocean. The input to 
this program is a data base, also described here, which 
contains the surface observations and coincident SMMR data. 
The evaluation of meteorological parameters using SMMR data 
is done by a user-supplied subroutine. Instructions are 
given for executing the program and writing the subroutine. 
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I. INTRODUCTION 


This report describes a program (PARAM) and data base that 
can be used for testing algorithms which are used to infer the 
values of meteorological parameters from data obtained by the 
SMMR instrument on NIMBUS-7. Each record in the data base 
combines a surface observation with coincident SMMR brightness 
temperatures. A coincidence is defined here as an observation 
by SMMR of the position where a surface observation was made 
where the time difference between the observations is less than 
1.5 hours. All surface observations referred to in this report 
were taken over the ocean. 

PARAM calls the user-supplied subroutine PARMS which may 
evaluate up to seven meteorological parameters from SMMR 
brightness temperatures. These parameters are freezing level, 
rain rate, wind speed, sea surface temperature, the amount of 
non-precipitating liquid water, the fraction of the field of 
view (FOV) of SMMR covered by rain, and the water vapor 
column density. Surface observations of wind speed, sea surface 
temperature, and water vapor column density are currently 
available in the data base. 

Sec. II describes the data base in more detail, and Sec. 

Ill tells how to use PARAM. Appendix 4 summarizes PARAM and 
all programs used to produce the data base. Listings of these 
programs are given in Appendices 5,6, and 7. 
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II. THE DATA BASF. 


The input to PARAM is a data base which consists of infor- 
station from surface observations and brightness temperatures 
obtained by SMMR. The surface observations come from tapes 
containing ship, buoy, and radiosonde reports (hereafter 
referred to as "station reports") which have been supplied 
by the National Climatic Center in Asheville, North Carolina. 
The brightness temperatures are read from SMMR CELL tapes. 

Station reports are included for which there exists 
a SMMR observation of the station's position at a time less 
than 1.5 hours before or after the time of the station 
report. The time limit of 1.5 hours is imposed in order to 
assure the validity of the comparison between the SMMR and 
station observations; this also prevents the same SMMR 
observation from being matched with more than one station 
report. 

There are other criteria for rejecting coincidences of 
station and SMMR reports. In the case of the ship and buoy 
reports, rejection occurs when both the wind speed and the 
sea surface temperature values are unusable. Radiosonde 
reports are rejected when the integration to determine the 
amount of water vapor above a point on the earth's surface 
fails. 

Brightness temperatures are available for the station's 
position and for the surrounding area as well. Each block 
on a CELL tape covers a box 780 km. on a side on the earth's 
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surface. For each coincidence between station and SMMR obser- 
vations, all of th brightness temperatures from the coincident 
CELL tape block are written into the database record along 
with the coincident station observation. 

Each database tape covers one month beginning with 
November 1978. The records on each tape are sorted in order 
of increasing SMMR observation time. 
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III. HOW TO USE PARAM 

PARAM is a main program, controlled by the namelist 
CNTL, which reads the database tapes, calls a user-supplied 
subroutine PARMS to evaluate meteorological parameters from 
the SMMR data, and produces printed output on request. The 
namelist CNTL, which is read in by PARAM on unit 5, is 
described in Appendix 1. Note that a user may restrict his 
analysis to a limited area on the earth's surface by specifying 
values for the namelist variables ALATL, ALATU, ALONL, and 
ALONU. 

The subroutine PARMS may evaluate up to seven meteoro- 
logical parameters: freezing level, rain rate, wind speed, 
sea surface temperature, the amount of non-precipitating 
liquid water, the fraction of the field of view (FOV) covered 
by rain, and the amount of water vapor. The parameters 
must always be passed back to PARAM in this order (see Appendix 
2). Appendix 2 also describes the way the brightness temper- 
atures are passed to PARMS. 

On return from PARMS, PARAM tests the value of the error 
flag, IERR. If I ERR is zero, an error has occurred in the 
calculation of the parameters, and PARAM will ignore the 
\ current record and read the next database tape record. If 

PARMS returns normally, IERR is set to a non-zero number. 

There are three parameters for which the values computed by 
PARMS can be compared with surface truth: wind speed, sea 
surface temperature, and water vapor. Water vapor is the only 
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parameter available from the radiosonde reports. Wind speed 
and sea surface temperature, but not water vapor, are given 
by the ship and buoy reports. 

After processing of the database tape has been completed, 
PARAM will print the average difference between SMMR values 
(as given by PARMS) and surface truth for any of the three 
parameters mentioned in the preceding paragraph (if the user 
has specified that they are to be evaluated) . The standard 
deviation of the difference and the number of usable measure- 
ments will be printed at the same time. Unusable measurements 
are those for which ground truth or SMMR parameters are not 
available. 

The user can make several different analyses with the 
same database tape by varying the namelist input. When PARAM 
finishes processing the tape, it checks for new namelist 
input. If more input exists, it rewinds the tape and process- 
ing starts from the beginning. 

Sample output from PARAM and the namelist input used to 
produce it are shown in Appendix 3. For each record of the 
database tape, there is a maximum of two lines of output. 

When two lines are printed, the first line refers to the first 
set of headings, and the second line belongs with the second 
set of headings. As described in Appendix 1, the user can 
specify that no printed output is to be produced, but that 
will not affect the printing of the average difference and 
standard deviation of the difference of the appropriate 
parameters. 
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APPENDIX 1. NAMELIST INPUT 


KOBS, KPARM, and KGRID must be specified by the user. Default 
values of the other namelist variables are given here. 


Variable 


KOBS 


KPARM 


Type Comments 

1*4 KOBS controls the type of station reports 

which are accepted. If KOBS is divisible 
by 2, ship reports are accepted. If it 
is divisible by 3, buoy reports are 
accepted. If it is divisible by 5, 
radiosonde reports are accepted. For 
example, if KOBS equals 10, only buoy 
reports are rejected. The minimum value 
of KOBS is 2 and the maximum is 30. 

1*4 KPARM tells PARAM which parameters are to 

be computed, so that it can set up the 
format to print them out. If KPARM is 
divisible by 2: freezing level is computed 
3: rain rate is computed 
5: wind speed is computed 
7: sea surface temperature is 
computed 

11: the amount of non-precipi- 
tating liquid water is 
computed 

13: the fraction of the field 
of view covered by rain 
is computed 
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17 : the amount of water 


vapor is computed 

The minimum value of KPARM is 2 and 
the maximum is 510,510. 

KGRID 1*4 KGRID * 1. Take brightness temperature 

from 156 km. cells. 

2. Take brightness temperature 
from 97.5 km. cells. 

3. Take brightness temperature 
from 60 km. cells. 

4. Take brightness temperature 
from 30 km. cells. 

LPRINT 1*4 Controls the amount of printed output 

LPRINT = 0: No printed output 
Is Print parameters 
2: Print parameters and station 
position and observation time. 
3: Print station position and 
observation time. 

The default is 2. 

The following four variables allow the user to reject reports 
from stations which are not within a specified area of the 
earth's surface. Note that the range of latitude is [-90 , + 90 ], 
and the range of longitude is [-180 , + 180*3 where longitude is 
positive to the east . 
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ALATL 

R # 8 

Lower 

limit 

of 

latitude. 



-90*. 




ALATU 

R*8 

Upper 

limit 

of 

latitude. 



+90*. 




ALONL 

R*8 

Lower 

limit 

of 

longitude 



-180*. 




ALONU 

R*8 

Upper 

limit 

of 

longitude 



+180*. 





Default Is 


Default is 


Default is 


Default is 


APPENDIX 2. HOW TO USE THE SUBROUTINE PARMS 


Brightness temperatures are passed to PARMS in the array 
TB which has 10 elements (five channels, two polarizations) . 
Each element occupies a halfword of storage, i.e. TB is an 
INTEGER* 2 array. The number of elements of TB which are 
filled depends on the grid size chosen: 


Number of 


KG RID Cell 

Size 

(km.) Brightness Temperature Frequency 

(GHz) 

1 

156 

10 

6.6, 10.7, 

18, 




21, 37 


2 

97.5 

8 

10.7, 18, 

21, 37 

3 

60 

6 

18, 21, 37 


4 

30 

2 

37 


The elements 

of TB 

are filled in the following 

order : 


Element 


Frequency (GHz) 

Polarization 


1 


6.6 

V 


2 


6.6 

H 


3 


10.7 

V 


4 


10.7 

H 


5 


18 

V 


6 


18 

H 


7 


21 

V 


8 


21 

H 


9 


37 

V 


10 


37 

H 



The brightness temperatures passed to PARMS are taken from the 
cell which covers the station's position. 


9 


The parameters calculated by PARMS are passed to PARAM in 
the array PARM. PARM is a REAL* 4 array with seven elements which 
are filled in the following order: 


Element 

1 

2 

3 

4 

5 

6 
7 


Parameter 

Freezing level (km) 

Rain rate (mm/hr) 

Wind speed (m/sec) 

Sea surface temperature (K) 

Non-precipitating liquid water (cm) 

Fraction of FOV covered by rain 

2 

Water vapor (gm/cm ) 


If any errors are encountered which render the results 
of the calculations in PARMS unusable, set the error flag 
IERR to zero. Conversely, if everything is normal, make sure 
that IERR is not zero. 

The first three lines of code in PARMS must be the following: 


SUBROUTINE PARMS (TB, PARM , IERR) 
DIMENSION PARM (7) 

INTEGER* 2 TB(10) 
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APPENDIX 3. SAMPLE OUTPUT FROM PARAM AND 
THE NAMELIST INPUT THAT GENERATED IT 
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APPENDIX 4. SUMMARY OF COMPUTER PROGRAMS 


A. Processing the Station Reports 

The first step in producing the data base tapes is the 
processing of the station report tapes from the National 
Climatic Center. A block diagram of the tape processing 
system is given in Figs. A4-1, A4-2, and A4-3. More detailed 
information can be obtained by consulting the listing in 
Appendix 5. 

If ship and buoy reports are being processed, MAIN calls 
BUSHIP. MAIN calls RADSON to process the radiosonde reports. 

Two subroutines are necessary because the radiosonde tape 
format is not the same as the ship and buoy tape format. 

BUSHIP calls the system subroutines POSN and FREAD to 
position the tape to the appropriate file and to read the tape. 
UNPKSB unpacks data concerning the station's position and the 
time it made an observation. These data are passed to OTIME 
which determines when SMMR could have observed the station's 
position in an interval of 48 hours centered about the station's 
observation time. OTIME assumes that SMMR is always turned on 
(see pt. B) . If SMMR observed the station within 1.5 hours of 
the station's observation time, WRITSB unpacks the meterological 
data on the statv.n report tape and writes out a record into 
the output file. 

This system possesses a restart capability. After every 
300 records are read from the station report tape, LTIME is 
called to check the amount of CPU and I/O time remaining. It 
does this by calling the system subroutine REMTIM. If less 
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BUSHIP 























YORB XORE 1 

















than 10 seconds of CPU or I/O time remain, the number of records 
read is printed and the job step is terminated. By setting the 
namelist variable ISKPRC to this number, the records which have 
been processed can be skipped the next time the station report 
tape is read. 

The block diagram for RADSON is similar to that for BUSHIP. 
The function NLEVEL determines the number of pressure levels 
at which the radiosonde made observations. 

After the tape is processed, the output file is sorted so 
that the records appear in order of increasing SMMR observation 
time. This is done to facilitate reading the SMMR CELL tapes. 
Processing of the station report tapes is completed by merging 
the sorted files into a single file. 

JCL used to run this system is given in Appendix 5 along 
with the listing. 

B. SMMR CELL Tape Processing System 

The purpose of this system is to find the block on a CELL 
tape which corresponds to the SMMR observation of the surface 
station. The output file created by this system is a database 
tape which is used as input to PARAM. Each record in this 
file consists of a station report and brightness temperatures 
from the corresponding CELL tape block. A block diagram of 
the system is shown in Figs. A4-4 and A4-5. 

MAIN calls NUFILE which positions the CELL tape to file 2, 
the first data file. GETSKP reads a station report, and GETCEL 
reads a block on the CELL tape. If the time of the SMMR 
observation, as given by the station report, is not within about 
2 minutes of the time corresponding to the CELL tape block, GETCEL 
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SMMR CELL Tape Processing System (continued) 


DISTNC 








will read forward in either the station or the CELL tape file 
until this condition is met. 

MAIN then calls INBLK to determine whether the area of 
the earth's surface covered by the CELL tape block includes 
the station's position. If the station is too far from the 
subsatellite track to be included in a CELL block, BACELL 
sets ICODE to zero, and the procedure is started again from 
the beginning, i.e. a new station report is read. If the 
station is not too far from the subsatellite track, INBLK 
performs a test to determine whether the SMMR observation 
of the station occurred in the current block or the next 
block or the preceeding block. If INBLK tells MAIN to read 
the preceeding block, BACKUP is called. If INBLK tells MAIN 
to read the next block, GETCEL is called. Once the new 
CELL block is read, the procedure described in this paragraph 
is repeated until INBLK decides that the current block covers 
the station's position. However, if INBLK decides to read 
forward on the CELL tape when it has previously decided to 
back up or vice-versa, MAIN will ignore the current ship record 
and read a new one in order to avoid getting into an infinite 
loop. The situation just described will occur when a block 
corresponding to the SMMR observation time as given by OTIME 
(see pt. A) is not written on the CELL tape either because 
SMMR was not on at that time or for some other reason. 

When INBLK decides that the appropriate block has been 
found MAIN calls RITEP to write a record onto the database as 
described in the first paragraph of this section. A listing 
of this system is given in Appendix 6. 
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C. PARAM 


The first task performed by PARAM (see listing in Appendix 7) 
is to check the namelist input. If the user has specified that 
the parameters calculated by PARMS are to be printed, PARAM sets 
up an execution-time format which allows for the printing of 
only those parameters and the corresponding surface truth (if 
any exists) . 

The next block of code begins with the statement which 
reads a data base record. If the ship is within the area 
specified by the user (see Sec. Ill and App. 1), the brightness 
temperatures are placed in the array TB as described in 
Appendix 2. Then PARAM calls the subroutine PARMS, and opera- 
tion of PARAM proceeds as described in Sec. III. 
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APPENDIX 5. LISTING OP STATION REPORT TAPE PROCESSING SYSTEM 


// EXEC FORTRANG 

//SOURCE ♦SYSIN DD DSN=ZMMHT . BUOY ♦ FORT > DISP*SHR 
// EXEC FORTRANH 

//SOURCE. SYSIN DD DSN=ZMMHT .BUOY. H. FORT »DISP*SHR 
// EXEC LINKGO, REGION. G0=148Nf0UT*X>TERM0UT=X 
//G0.FT05F001 DD * - 

SINPUT KIND*1*ISKPRC*0, MAXREC=1000000>LPRINT*F»ATAPE*'NCC001 ' «• SEND 
//GO.FTOSFOOl DD SYSOUT*X — • 

//GO ♦ F-T09F00 1 DD DSN*ZMhHT . N0V78DAT r DISP=SHR 

//G0.FT15F001 DD LABEL*< 1 ,NL , » IN ) , - * 

// DCB* ( LRECL*1 40 1 BLKSIZE*1 400 1 RECFM=FB t DEN=3 ) * 

// VOL*SER=BUOYlf UNIT* <9TRACKr» DEFER) 

//G0.FT17F001 DD UNIT*DISK. DSN=ZMMHT . N0U781 1 .US.DATA> 

// DISP*(NEUf CATLG) » SPACE* (CYL» 10) » 

// DCB* < RECFM=VBS »LRECL=92 * BLKSIZE=924 ) 

//G0.FT11F001 DD DSN=ZMMHT.ELTS.N0U78.DISP=SHR . __ 


//SORT EXEC' S0RTD»PARM*'C0RE=250K»0PT=D / 

//SORTIN DD DSN=ZMMHT.N0U7811 .US. DATA f DISP* (OLD f DELETE) 
/ /S0RTWK01 DD UNIT=SYSDA»SPACE=(CYL* 2) 

//S0RTUK02 DD UNIT*SYSDA r SPACE* (CYL»2) 


/ /S0RTUK03 DD UNIT*SYSDA» SPACE=(*CYL »2) 
/ /S0RTWK04 DD UNIT=SYSDArSPACE=(CYL»2) 


3 .7 ."i 


/ /S0RTWK05 DD UNIT=SYSDA,SPACE=<CYL»2) 

//SORTOUT DD UNIT=DISK » DSN=ZMMHT .N0V7811 .SORTED.DATA? 

// DISP*<NEW»CATLG) fSPACE=(CYLf 10) » • r 

// DCB* ( RECFM=VBS » LRECL=92 r BLKSIZE=924 ) r . 

SORT FIELDS=<37>4f BIr A»41 »4f BI f A. 4 B if 8rBI r A) r SIZE*E1000009 
RECORD TYPE*V» LENGTH* (92f92*92) 


END 


Fig. A5-1. JCL used to process a ship tape and sort the output. 
This is run as one job. 


Subroutines appear in alphabetical order after the main 
program. OTIME, GKEP, LONG IT. GETTIM, XORB, and YORB are listed 
in the report "Finding the Times that SMMR Observed a Ship" 
and do not appear here. 
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//MERGE EXEC S0RTDjPARM«'C0RE-250K>0PT-D' 

//S0RTIN01 DD UNIT-DISK* DSN-ZMMHT ♦ NQU78 1 1 ♦ SORTED ♦ DAT A r DISP-OLD 
//S0RTIN02 DD UNIT-DISK » DSN-ZMMHT . N0U7821 . SORTED . DATA , DISP-OLD 
//S0RTIN03 DD UNIT-DISK t DSN-ZMMHT ♦ N0U7831 . SORTED . DATA t DISP-OLD 
//SORTOUT DD UNIT-(9TRACK» » DEFER) »LABEL-( 1 »NL» *OUT> »DISP=<NEW.KEEP> 
// DCB-(RECFM*UBS f LRECL-92 1 BLKSIZE-924 1 DEN-3 ) t 

// VOL-SER-MHT 004 

//SYSIN DD * ^ 

MERGE FIELDS»<37,4*BI*A*41*4*BI*A*45»8»BI*A) , SIZE-E1000009 
RECORD TYPE-V t LENGTH- (92*92*92) 

end 


Ei<J« A5-2. JCL used to merge sorted ship, buoy, and radiosonde 
files. 
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APPENDIX 6. LISTING OF SMMR CELL TAPE PROCESSING SYSTEM 


//STEP1 EXEC FORTRANG >OUT*A 

//SOURCE. SYSIN DD DSN*ZMhHT .LHB .FORT t DISP*SHR 
// DD DSN*ZMMHT .CELL. FORT >DISP=SHR 

//STEP2 EXEC LOADER »REGI0N=19SK >PARM* ' SIZE=150K ' »OUT=A t TERMOUT=A 
//£T06F00l DD SYSQ UT»A 

//FTllFOOl DD UNIT*( 9TRACK» > DEFER) * ' ~ 

// DCB«(LRECL»15120> BLKSIZE-15120 »RECFM=FB »DEN=3 ) > 

// LABEL*(2»NLr.» IN) t U0L*SER*L8329, 

// D ISP»< OLD > KEEP) 

//GO .FT14F001 DD UNIT»(9TRACK» * DEFER) r LABEL* < 1 »NL» fOUT ) t 
// DCB= < RECFMr'VBS » LRECL*6496 > BLKSIZE*6500 » DEN=3 ) » 

// DISP“( NEW* KEEP) » V0L*SER=MHT001 

//GO .FT20F001 DD UNIT* < 9TRACK > » DEFER) f LABEL* < 1<-NL» »IN) r 
// DCB* ( RECFM=UBS * LRECL*92 » BLKSIZE=924 » DEN=3 ) » 

// DISP«<OLD* KEEP) r U0L=SER=MHT004 . . .... 


Fig. A6-1. JCL used to run SMMR CELL tape processing system. 
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APPENDIX 7. LISTING OF PARAM 


// EXEC FORfGl 

//SOl'RCE ♦ SYSIN DO DSN-ZMMHT.LIB2.F0RT<PARAM) rDISP-SHR 
// EXEC L INKGO t REG I ON . G0»98K ? OUT *A r TERMOUT «A 
//G0.FT06F001 DD SYSOUT-A 

//GO*FT11F001 DD . UNIT*(9TRACK» » DEFER) » LABEL* ( 1 »NLf » IN) » 
// DISP* ( OLD f KEEP) » V0L*SER“MHT001 r 

// DCB-<RECFM«VBS>LRECL«A49A>BLKSIZE»6500> DEN-3) 

//G0.DATA5 DD * 

ICNTL KOBS*2» KPARM*510510> KGRID*1 t SEND 


Fig. A7-1. JCL used to run PARAM. 
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